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We develop a full microscopic replica field theory of the dynamical 
transition in glasses. By studying the soft modes that appear at the 
dynamical temperature we obtain an effective theory for the critical 
fluctuations. This analysis leads to several results: we give expres- 
sions for the mean field critical exponents, and we study analytically 
the critical behavior of a set of four-points correlation functions from 
which we can extract the dynamical correlation length. Finally, we 
can obtain a Ginzburg criterion that states the range of validity of 
our analysis. We compute all these quantities within the Hypernet- 
ted Chain Approximation (HNC) for the Gibbs free energy and we 
find results that are consistent with numerical simulations. 

glass transition | mean-field theory | dynamical heterogeneities 

Introduction. Dynamical heterogeneities in structural glasses 
have been the object of intensive investigations in the last 
15 years [1]. The early Adams-Gibbs theory of glass forma- 
tion was based on the concept of cooperatively rearranging 
regions, whose size becomes larger and larger when the glass 
region is approached. Such large cooperatively rearranging 
regions imply the existence of dynamical heterogeneities char- 
acterized by a large correlation length. Large scale dynamical 
heterogeneities are expected to be present in any framework 
where glassiness is due to collective effects: they are indeed 
the smoking guns for these effects [2, 3, 4, 1]. Therefore, it is 
not a surprise that two popular approaches to glasses, Mode 
Coupling Theory (MCT) [5] and the replica method [6, 7], 
both agree with the Adams-Gibbs scenario and predict large 
scale dynamical heterogeneities with a dynamical correlation 
length that diverges at the transition to the glass phase. This 
qualitative prediction is very interesting, but in order to make 
further progresses it would be important to get quantitative 
predictions, that can be compared with numerical simulations 
and with experiments. 

At the mean field level, where both thermodynamic and 
dynamic aspects can be solved exactly, it is found that the 
replica and Mode-Coupling approaches are intimately related. 
The study of spherical p-spin models, where dynamics is ex- 
actly described by a schematic MCT equation and equilibrium 
display glassy phenomena related to Replica Symmetry Break- 
ing (RSB), shows how the glass transition described by MCT 
is related to the emergence of metastable states in equilib- 
rium [8, 9]. That basic observation, made more then 20 years 
ago by Kirkpatrick, Thirumalai and Wolynes [10] , opened the 
way to the application of the mean field theory of spin glasses 
to the physics of supercooled liquids and glasses [11, 12, 13]. 
Despite this clear relation at the level of mean field schematic 
models, when one tries to apply the mean field theory to real- 
istic models of simple liquids [14, 15, 6, 5, 7] approximations 
are mandatory, and because of that the connection between 
statics and dynamics becomes more difficult to establish. It 
has been shown by Szamel [16] that under suitable approx- 
imations, similar to the one of MCT, the long time limit of 
the MCT equations could be derived from a replicated liq- 
uid theory. Unfortunately this leads to expressions that are 
not variational and one cannot get an approximation for the 
free energy from the computation. Using instead standard liq- 



uid theory approximations within replica theory [14, 15, 6, 7], 
one finds strong discrepancies between predictions from MCT 
and replicas which become particularly pronounced in large 
dimensions [7, 17]. 

Besides this consistency problem, in finite dimensions one 
would like to compute the corrections due to fluctuations 
around the mean field approximation. When this program 
is carried out, one finds that there are two important sources 
of corrections to the mean field scenario. The first correc- 
tions originate from critical fluctuations that become impor- 
tant around the glass transition below the upper critical di- 
mension, as in any standard critical phenomenon [18, 19]. The 
second corrections are non-perturbative phenomena related to 
activated processes. They can be taken into account by a phe- 
nomenological approach, leading to a number of predictions 
that are in good agreement with experiment [11]; however, 
the theoretical foundations of this approach are still contro- 
versial [20] and alternative (but possibly related) phenomeno- 
logical descriptions of activated relaxation in glasses have been 
developed, mostly based on the concept of dynamical facilita- 
tion [21]. 

In this paper we will only consider critical fluctuations 
around mean field, so we will not take into account acti- 
vated processes. Critical fluctuations have been previously 
described within MCT [22, 18, 23, 24]. However, field the- 
oretical methods are not yet under complete control in the 
context of dynamics, and it is therefore extremely important 
to set up a static replica field theoretical description of dy- 
namical heterogeneities, in such a way that well-established 
equilibrium field theory methods such as the renormalization 
group can be applied to the glass transition problem. This is 
what we achieve in this paper. We obtain a low-energy effec- 
tive action that describes critical fluctuations on approaching 
the glass transition, whose coupling constants are obtained di- 
rectly from the inter-particle interaction potential using stan- 
dard liquid theory. This allows us to compute prefactors to the 
singular behavior of physical observables in the mean field ap- 
proximation, such as the correlation length or the four-point 
correlation functions. In addition, we show that an impor- 
tant characterization of dynamics, the MCT exponents, can 
be obtained within the static replica framework. Using the 
well established HNC approximation of liquid theory, we per- 
form explicit computations for Hard and Soft Sphere models 
and Lennard- Jones potentials and we obtain good agreement 
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with available numerical data. Finally we introduce a quanti- 
tative Ginzburg criterion defining a region where perturbative 
corrections to mean field theory can be neglected. 

Dynamical heterogeneities. In the following we consider a sys- 
tem of N particles in a volume V interacting through a pair- 
wise potential v(r) in a D dimensional space. The dynamical 
glass transition is characterized by an (apparent) divergence of 
the relaxation time of density fluctuations, that become frozen 

in the glass phase. If p(x,t) = YLi=i $(. x ~ Xi W) * s ^ ne l° ca l 
density at point x and time t and p = {p(x, t)) its equilibrium 
average, the transition can be conveniently characterized us- 
ing correlation functions. Consider the density profiles at time 
zero and at time t, respectively given by p(x,0) and p(x,t). 
We can define a local similarity measure of these configura- 
tions as 

6{r,t) = J dxf(x)p(r + ^,t)p(r~^,0)~p 2 , [1] 

where f(x) is an arbitrary "smoothing" function of the den- 
sity field with some short range A. In experiments, f(x) could 
describe the resolution of the detection system and can be for 
instance a Gaussian of width A. 

Let us call C(t) = V' 1 J dr{C(r,t)) the spatially and 
thermally averaged correlation function. Typically, on ap- 
proaching the dynamical glass transition Td, C(t) displays a 
two-steps relaxation, with a fast "/3-relaxation" occurring on 
shorter times down to a "plateau", and a much slower "a- 
relaxation" from the plateau to zero [5]. Close to the plateau 
at C(t) = Cd, one has C(t) ~ Cd + At~ a in the /3-regime. 
The departure from the plateau (beginning of a-relaxation) 
is described by C(t) ~ Cd — Bt b . One can define the a- 
relaxation time by C(r a ) = C(0)/e. It displays an apparent 
power-law divergence at the transition, r Q ~ \T — Td| -7 . All 
these behaviors are predicted by MCT [5]. In low dimensions, 
a rapid crossover to a different regime dominated by activa- 
tion is observed and the divergence at Td is avoided; however, 
the power-law regime is the more robust the higher the dimen- 
sion [25, 26] or the longer the range of the interaction [27]. 

It is now well established, both theoretically and experi- 
mentally, that the dynamical slowing is accompanied by grow- 
ing heterogeneity of the local relaxation, in the sense that the 
local correlations C(r, t) display increasingly correlated fluc- 
tuations when Td is approached [2, 28, 3, 1]. This can be 
quantified by introducing the correlation function of C(r,t), 
i.e. a four-point dynamical correlation 



G 4 (r,i) = <C(r,t)C(0,t)> - <C(r,t))<C(0,*)> 



[2] 



The latter decays as Gi(r,t) ~ exp(— r/£(t)) with a "dynam- 
ical correlation length" that grows at the end of the /3-regime 
and has a maximum £ = £(i ~ r a ) that also (apparently) 
diverges as a power-law when Td is approached. 

MCT [5] and its extensions [22, 18, 23, 24, 29, 30] give 
precise predictions for the critical exponents. However, as dis- 
cussed in the Introduction, this dynamical transition can be 
also described, at the mean field level, in a static framework. 
This has the advantage that calculations are simplified so that 
the theory can be pushed much forward, in particular by con- 
structing a reduced field theory and setting up a systematic 
loop expansion that allows to obtain detailed predictions for 
the upper critical dimension and the critical exponents [19]. 
Moreover, very accurate approximations for the static free en- 
ergy of liquids have been constructed [31], and one can make 
use of them to obtain quantitative predictions for the physical 
observables. This is the aim of the rest of this paper. 



Connection between replicas and dynamics. In the mean field 
scenario, the dynamical transition of MCT is related to the 
emergence of a large number of metastable states in which the 
system remains trapped for an infinite time. At long times in 
the glass phase, the system is able to decorrelate within one 
metastable state. Hence we can write 



<<7(r,t-K»)) = / dxf(x)(p(r + 



pir 



where (•) m denotes an average in a metastable state, and the 
overline denotes an average over the metastable states with 
equilibrium weights. 

The dynamical transition can be described in a static 
framework by introducing a replicated version of the sys- 
tem [32, 14]: for every particle we introduce m— 1 additional 
particles identical to the first one. In this way we obtain m 
copies of the original system, labeled by a = 1, . . . ,m. The 
interaction potential between two particles belonging to repli- 
cas a, b is v a b(r). We set v aa (r) = v(r), the original potential, 
and we fix v ab (r) for a 7^ b to be an attractive potential that 
constrains the replicas to be in the same metastable state. Let 
us now define our basic fields that describe the one and two 
point density functions 



[4] 



Pab( X >V) = Pa(x)pb(v) ~ pa(x)S ab S(x - y) 



To detect the dynamical transition one has to study the 
two point correlation functions when v a b(r) — > for a 7^ 
b, and in the limit m — > 1 which reproduces the original 
model [32, 14]. In this limit, the two-replica correlation func- 
tion is, for a 7^ b: 

{C ab (r)) = J dxf(x)(p a (r+^)p b (r-^))-p 2 . [5] 

Because of the limit V ab (r) —> 0, the two replicas fall in the 
same state but are otherwise uncorrelated inside the state, 
therefore we obtain {C a b{r)) = (C(r,t — > 00)) which provides 
the crucial identification between replicas and dynamics. Sim- 
ilar mappings can be obtained for four-point correlations. 

Replica field theory for the dynamical transition. We intro- 
duce for convenience an external field v a {x) (that derives from 
a space-dependent chemical potential), in such a way that the 
density correlation functions can be obtained by taking the 
derivative of the free-energy with respect to it [31]. The free 
energy is defined as the logarithm of the partition function, 
and its double Legendre transform defines the Gibbs free en- 
ergyr[{p a (x)},{p$(z,j/)}] [31,33]: 



[ Pab( X yy) 

Pab ( x ,y)i n ' 



Pab( x >y) + p*( x )pb(y)] + / dx P*( x ) [ ln Pa(aO - 1] 

a •> 

{-If 



2n 



+ E 

X ' ' ■ Pa„(x n )h a „ ai (x n ,Xl) +r 2 PI 



dXl ■ ■ -dXn p ai (xi)h ai a 2 (Xl,X2)X 



where h ab (x,y) = pf b '(x,y)/p a (x)p h (y) - 1 and r 2P i is the 
sum of 2-line irreducible diagrams [33]. The average values 
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of the fields in Eq. [4], namely p a (x) and p ab (x,y), can be 
obtained by solving the saddle point equations 



ST[{pa},{p%} 



Po.b( x -V) 



[7] 



and similarly for p a (x). Here we consider a homogeneous liq- 
uid, hence pa(x) = p. 

We have to assume at this point that a mean field approx- 
imation of the free energy is available, that we shall use as the 
starting point of our computations. Within this approxima- 
tion, we want to study the behavior of~p a7ib (x,y) in the double 
limit m — > 1 and v a ^b — > 0, which signals the dynamical tran- 
sition: if T > T d , then p a ^ b (x,y) = p 2 while if T < T d a non 
trivial off-diagonal solution persists in the limit v a ^b — > 0. At 
the mean field level, the appearance of the non trivial solution 
is a bifurcation phenomenon so that, if we come from below 
the transition and we define e = Td — T, we have for e — > 0: 

Pa#b( 2; .y; e ) = P 2 9( x ~ v) + 2 P 2 Ve « ko{x-y) , [8] 

where ko(x) is normalized as J dxko(x) 2 = 1, and k is a 
constant. From the saddle point equations [7] we obtain 
that the Hessian matrix for the off-diagonal elements, i.e. for 
a / 6, c / d 



M ab . cd (x 1 ,X 2 ;X 3 ,X i ) = 



* 2 r[{M,{pff}] 



considered as a kernel operator both in standard and replica 
space develops a zero mode at Td . This means that if the tran- 
sition is approached from below, the fundamental eigenvalue 
of this operator is proportional to i/e due to the bifurcation- 
like phenomenology. Moreover the eigenvector corresponding 
to it is ko(x — y). 

Exploiting the replica symmetry of the saddle point solu- 
tion Eq. [7], the most general form of the Hessian matrix is 
given by 

\ 8 ac 8bd + 5adSbc\ 

M ab;C d{xi,X2;x3,X4) = Mi [ ] + 



+M 2 [' &"= + *<-* + &° + *m ] + M3 



[101 



where Mi, M 2 and M3, depend on xi, . . . ,£4. From this one 
can show that, because the zero mode ko(x — y) is indepen- 
dent of the replica indices, in the replica limit m — > 1 it is 
an eigenvector of the kernel operator Mi . To study the cor- 
relation functions for the fields in Eq. [4] we can produce 
a power series expansion of the Gibbs free energy in terms 
of the fluctuation of the field p a 2 ^ b (x, y) from its saddle point 

value. Defining the field Ap ab (x,y) = p%}(x,y) -p ab (x,y), 
we can expand the Gibbs free energy up to the third order. It 
is convenient to define pi and qi as the momenta conjugated 
to the half sum and the difference of the spatial arguments of 
Ap ab (xi, yi). Using translation invariance we write the replica 
action in Fourier space as 



r[{A Pab }] = r[{p ab }]- 

1 ^ f dpdq 1 dq 2 



&Pab(p, qi)M ( a p b ) d (qi, q 2 )Ap cd (~p, q 2 ) + 



2 ^ J (2tt) 3D 



1 f dpdp'dqidq 2 dq 3 r , , . 

6 1^ / 7^5 L ab , cd , ef (p,p;q 1 ,q 2 ,q i )x 

x Ap ab (p, 5i)Ap cd (p', q 2 )Ap ef (-p -p',q 3 ) 



Because of the zero mode of the Hessian matrix, the connected 
correlation function of Ap ab (x,y) shows critical fluctuations 
at the transition. 

To make the connection with the dynamical correlation, 
we define an overlap function among replicas, q a b(r) as in 
Eq. [1] substituting the configurations at time and t by 
replicas a and b. We expect that all the critical fluctuations 
of q a b(r) can be captured by a projection on the zero mode, 
leading from Eq. [11] to an effective action. We can study 
the fluctuations of q a b(r) for generic functions /, by perform- 
ing a Legendre transform of Eq. [11]. However the results 
are quite involved and here for clarity we will first consider 
the simplest case where f(x) = ko(x). Of course, this is not 
a practical choice for numerical simulations or experiments 
because ko is quite difficult to measure, however the theoret- 
ical computations are much simpler in this case. Later on we 
will show that any other choice of / leads to the same results 
for the critical quantities, and it only affects the prefactor of 
the correlation functions. The projection onto the zero mode 
can be done by choosing Ap a b(x,y) = ko(x — y)4> a b an d 
substituting this in Eq. [11]. The field (j> a b{x) is the com- 
ponent of the overlap along the zero mode, and we perform a 
perturbative expansion at small momentum p. The effective 
replica field theory that arises is equivalent to a Landau-like 
gradient expansion along the critical modes: 



+m 2 Y i 



^<t>ab{p) 



+ m 3 



a^b 



[12] 



+ 



Wo 



( ^I7r^4>ab(p)<t>ab{p')<i>ab(-P - P) , 
J a±b [ - Z7V> 



Eq. [12] is the effective low-energy replica field theory we 
will use to compute the critical properties of the system. All 
its coefficients can in principle be computed from the micro- 
scopic details of the systems once an approximation for F is 
available. In fact they can be given explicit expressions as 
functions of derivatives of the Gibbs free energy and of the 
zero mode, which both derive from the interaction potential 
(see Supporting Information). 

Correlation functions, correlation length and critical expo- 
nents. The effective replica field theory in Eq. [12] can be 

used to compute the MCT parameter A. This quantity is re- 
lated to the MCT critical exponents that control the approach 
to the plateau by the relation 



A 



rd 



_J_ _ r(i + b) 2 

T(l - 2a) ~ T(l + 26) 



[131 



In addition, the exponent that controls the growth of the re- 
laxation time r a ~ |T-T d |~ 7 is given by 7 = 1/ (2a) + 1/(26). 
Although A is a dynamical parameter, it has been explicitly 
shown recently in disordered mean field models and it can be 
argued on general ground [34, 35], that this parameter can be 
related to a ratio of 6 point static correlation functions com- 
putable in the replica field theory that we have just derived. 
In this scheme the exponent parameter is given by 



[11] 



Wi 



[141 
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Moreover, the field theory above can be used at the Gaus- 
sian level in order to obtain the correlation functions of the 
overlap. The analysis of the quadratic part of Eq. [12] shows 
that the correlation length is controlled by the diagonal part, 
being and mz finite at the transition. The result is 



-1/4 



ft 



[151 



and it corresponds to the divergence of the dynamical corre- 
lation length £(t) in the /3-regime [22, 23, 24]. 

Moreover we can compute in details the critical behav- 
ior of many possible dynamical four-point functions, that are 
identified with different matrix elements of the inverse of the 
Hessian matrix in Eq. [9], see [19]. Here we give the re- 
sults for the simplest one, the so-called in-state, or thermal 
susceptibility, that is given by 



i (r,t) = -E [{C(r,t)C(0,t)) - {C(r,t)){C(Q,t)) 



[161 



where in the equation above Eo[-] has to be intended as the 
average over the initial positions of the particles, while (•} is 
an average over different trajectories (i.e. over the noise for 
Langevin dynamics, or over the initial velocities for Newton 
dynamics). In the long time limit, this quantity is one of the 
critical contributions to the Gi(r,t) in Eq. [2], and it can 
be computed directly from the replica field theory above [19] . 
Here we had to generalize the calculation of [19] to take into 
account the structure of the zero mode and the presence of 
the smoothing function f(x). The result is 



Gth(p) = 



G e 



-1/2 



i + £V 



G = I/ 
V J 



(271-)- 



r f(-q)k (q). [17] 



We obtain that the correlation length and its prefactor are 
not dependent on the function f(x) and always given by 
Eq. [15]. The only dependence on f(x) of the four-point 
function is in the prefactor Go- The full four-point correla- 
tion [2] is known to display a doubled singularity with re- 
spect to [17]. In fact, with the choice f(x) = ko(x) one finds 
G 4 (p) = Gth(p) - (ma + m 3 )G th (p) 2 [19]. For generic f(x), 
the computation of the prefactor is more involved and will not 
be presented here. 

A Ginzburg Criterion. All the calculations above are based on 
the assumption that a mean field approximation of the free 
energy of the system is given. From this, we derive the ef- 
fective Landau field theory Eq. [12]. From its coefficients, 
we extracted all the mean field critical exponents, as well as 
microscopic expressions for the prefactors. Now we can check 
whether loop corrections to the effective field theory affect 
strongly the mean field predictions, by means of a Landau- 
Ginzburg computation. In other words we want to see whether 
the loop corrections to the bare correlation function are small. 
In principle we should take the field theory derived above, and 
then we should compute the first non trivial loop diagrams 
which give the first correction to the propagator in replica 
space. This computation is quite involved because we have to 
deal with replica indices. However it has been shown in [19] 
that the leading divergent behavior of the above field theory 
can be mapped to the one of a scalar field in a cubic potential 
with a random field 

5(0 = i / dxcb(x)(-aV 2 + p^+Sm(g,A))^x) 

or r [18] 

+ | / dxcj > 3 (x) + / dx(h (x) + 5h{g,A))cj>(x). 



where the random field has zero mean and correlation 
ho(x)ho(y) — A5(x — y), and the coupling constants are given 
by g = W2 — w\ and A = — 7712 — 7713. 

The terms Sm(g, A) and 5h(g, A) are counterterms needed 
to enforce that the critical point is not shifted by loop cor- 
rections. By computing the first one- loop diagrams and by 
imposing that the relative correction is small with respect to 
the bare quantity, we arrive to the following Landau-Ginzburg 
criterion 

-D 



1 > Gif 

where the (dimensional) Ginzburg number is given by 



Gi = 



4(4^)0/2 



19 



[201 



This computation is correct only below the upper critical di- 
mension D u = 8. For D > D u , the theory is divergent in 
the ultraviolet and the Ginzburg number depends on the mi- 
croscopic details, but the critical exponents coincide with the 
mean field ones. A similar calculation in the framework of 
MCT has been carried out by Szamel [36]. 

Results in the HNC approximation. Up to now the calcula- 
tions were very general and the results above hold for any 
given approximation of the replicated free energy functional 
that displays the correct mean field glassy phenomenology. 
One of the advantage of our static approach is indeed that it 
can be systematically improved by considering more accurate 
approximations of F. 

Here we report results obtained from the replicated HNC 
approach, that amounts to neglecting the F2pi term in 
Eq. [6], and has been shown to give the correct glassy phe- 
nomenology at the mean field level [14, 15]. Applying the 
formulae above, we find that in the HNC approximation the 
parameter A is given by 



A = 



■ fax 



2 ^/l2^0 3 to[l- 



P Ac(q)] 



[211 



where g(x) = p a ^ b (x) / p 2 , Ac(q) = c aa (q) - c a ^b{q), and the 
direct correlation function c a b(q) is related to h a b(q) by the 
replicated Ornstein-Zernicke relation [14] . Similar expressions 




Fig. 1. The zero mode ko(cf), the structure factor S(q) and the non-ergodicity 
factor f(q) for Hard Spheres at the dynamical transition = 1.176 in the HNC 
approximation. 
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can be obtained for all the other coefficients, see Supporting 
Information. 

To produce concrete numerical results we have solved nu- 
merically the HNC equations by standard methods [14] for 
a large variety of systems in D = 3. In particular we have 
considered 

• Hard Spheres (HS): v(r) — for r > ro and v(r) = oo 
otherwise. 

• Harmonic Spheres (HarmS): v(r) — e(ro — r) 2 6(ro — r). 

• Soft Spheres (SS-n): v(r) = e(r /r)", with n = 6,9, 12. 

• Lennard- Jones (LJ): v(r) = 4e [(r /r) 12 - (ro/r) ] 

• Weeks-Chandler- Andersen (WCA): 

v(r) = 4e [(r /r) 12 - (r /r) 6 + 1/4] 0(r o 2 1/6 - r) 

In all cases we fix units in such a way that ro = 1, e = 1 
and the Boltzmann constant fee = 1. For HS and SS, tem- 
perature is irrelevant (for SS the only relevant parameter is a 
combination of density and temperature, hence we fix T = 1 
for convenience), and we study the system as a function of 
density to determine the glass transition density pd- For the 
other systems, we studied the transition as a function of both 
density and temperature. 

In order to obtain numerically the zero mode we have used 
the definition in Eq. [8], and estimated it by the numerical 
derivative of g(r) with respect to *Je when e — > 0. A plot of 
the zero mode for HS is in Fig. 1. Interestingly we find that 
the zero mode has the same structure in Fourier space as the 
static structure factor S(q) and the non-ergodicity parameter 
f(q), which is the Fourier transform of the long time limit of 
Eq. [ 1 ] in the glass phase [5] . This finding offers a rationaliza- 
tion of the common practice of concentrating on momenta of 
the order of the peak of S(q) in the study of glassy relaxation. 

From the zero mode we can compute all the coefficients 
of the effective action from which we obtain the physical 
quantities. In particular, we can compute the prefactor £o 
of the growth of the correlation length and the Ginzburg 
number. Moreover, we have computed the prefactor Go of 
the in-state susceptibility Eq. [17] using a box function 
f(x) = (2A)- D/2 nf =1 9(A 2 -x 2 a ) where 9{x) is the Heaviside 
step function and A = O.lro. All the results are collected in 
tables 1 and 2. 

The value of A we find is almost the same for all inves- 
tigated systems but it is not consistent with the result of 
MCT [5] and with numerical results for these systems. Note, 
moreover, that the location of the critical point predicted by 
HNC is different from the one of MCT: e.g. for HS, HNC pre- 
dicts p d = 1.169 while MCT predicts p d = 0.978 [5]. This is 
an example of the fact, already mentioned in the introduction, 
that different approximation schemes lead to different results. 



Another example of this problem is obtained by comparing 
the results for LJ and WCA at p = 1.2, 1.4 (table 2) with 
MCT and numerical data reported in table 1 of Ref. [37]. 
The most interesting numerical result is the Ginzburg num- 
ber. We predict that (perturbative) corrections to mean field 
results in D — 3 should remain small as long as the dynamical 
correlation length is smaller than ~ 1. Note that a different 
Ginzburg criterion for the validity of MCT, based on a phe- 
nomenological approach, has been derived in [20]: the results 
of that analysis also suggest that corrections to mean field will 
appear when the correlation length is ~ 1. 

Unfortunately, not many data for the critical behavior of 
four-point correlations in the /3-regime are available [38, 39]. 
It would thus be very interesting to get high precision simu- 
lation data in the /3-regime. 

Conclusions. We have studied in details the replica field the- 
ory for the dynamical transition in glasses. By using the HNC 
approximation we have computed many physical observables 
directly from the microscopic expression of the interaction po- 
tential. First of all we provided a way to compute the Mode- 
Coupling exponent parameter A. The numerical values ob- 
tained are in good agreement with the experimental and nu- 
merical estimates. Moreover we have computed the prefactor 
of the correlation length at the transition, together with the 
prefactor of the in-state four-point correlation function. Fi- 
nally we have closed self-consistently our analysis by looking 
at the loop corrections to the mean field quantities in order to 
produce a Ginzburg criterion that states how close we have to 
be to the dynamical transition in order to see deviations from 
mean field theory. We found that the range currently acces- 
sible to numerical simulations in three dimensions is close to 
the point where such corrections should become important. 
Of course, non-perturbative corrections (activated processes) 
are not included in our analysis, but they are responsible for 
strong deviations from the MCT regime when the transition 
is approached. 

Our analysis is quite general because it relies only on the 
assumption that the approximation scheme used for the Gibbs 
free energy shows the correct mean field glassy phenomenol- 
ogy. Hence, it can in principle be repeated in different ap- 
proximation schemes in order to go beyond HNC and obtain 
more accurate expressions for physical quantities. 
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disregard the excited modes which are finite and take into ac- 
count only the projection of the dynamical field p' 2 ^ on the 
zero mode. Practically this is the same as putting to infinity 
the masses relative to the projections of the dynamical field 
on the excited states of the kernel operator Mi. By doing this 
we can produce a gradient expansion for the replicated Gibbs 
free energy. The simplest way to do this is to impose that 
the fluctuations of the dynamical field from the saddle point 
solution are proportional to the zero mode 



&Pab(x,y) = 



x + y 



ko(x - y) 



[221 



By doing this, the expressions for the coefficients of the ef- 
fective action for the critical fluctuations can be computed 
straightforwardly. Let us consider first the expression for a 
and \x. They come along in this way. The kernel operator Mi 
has a ground state eigenvalue Ao(p) = py/c + ap 2 + 0(p 4 ) . For 
small momentum (which means that we look at the correla- 
tion of two fluctuations of the dynamical field that are at a 
very large distance) the expressions for p and a can be com- 
puted using perturbation theory for the eigenvalue problem 
for the kernel Mi where the small perturbative parameter is 
exactly the momentum p. The final expressions are given by 



.. f d D gd D fc, . . d ar(p ), 
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[241 



where the zero mode is supposed to be normalized. In the 
same spirit the two other masses m;, i = 2,3, are given by 



rrii — lim 

£->0 



A D qA D k ) 
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Supplementary Information 

The Supporting Information text is organized in three parts. 
In the first one we give a sketch of the line of reasoning that 
leads to the effective action used to describe the dynamical 
transition and we give all the expressions for the coefficients 
of the same action in terms of the interparticle potential. First 
we do this in a generic framework without specifying the ap- 
proximation used to compute the Gibbs free energy and then 
we give the expressions in the HNC case. The second part 
of the present text is devoted to the Ginzburg criterion: wc 
describe the guidelines of the computation by showing which 
diagrams have been taken into account to produce the first 
correction to the bare four-point function that has been given 
in the main text. The last section contains some details on 
the numerical calculations and it is useful just to understand 
how our results can be improved numerically. 

Coefficients of the replica Gibbs free energy. As in the main 
text, we assume that the glassy phenomenology manifests it- 
self in the singular behavior of the off-diagonal field p ai tb{x, y) 
that has a diverging derivative with respect to temperature 
when the critical point is approached. This implies that the 
Hessian (or mass) kernel operator develops a zero mode. Ac- 
tually, we remember here that due to the replica symmetry of 
the saddle point we have that only one {i.e. Mi) of the three 
kernel operators Mi, M2, and M3 has a zero mode. This 
implies that the field p^ can be decomposed on the eigen- 
vectors of Mi. Because we want to give the expressions for 
the diverging part of the correlation function, we can simply 



At this point it is clear how the expressions for the two cubic 
coefficients wi and W2 can be obtained; defining 

->->ab;cd\ef [Xl^ - . . , Xq ) — t^t t^T t^T 

Sp a i! (xi, x 2 )8p y J (x 3 ,x 4 )Spl/ (x 5 , x 6 ) 

[26] 

then they are given by the following expressions 

101,2 = J d D xi, ... d D x 6 ko{x! — X2) ... k (x 5 — xe)Wi 7 2 

[27] 

where 

Wl = L a b,bc,ca — 3L a b, a c,bd + 3L aCi (,c,de — L a b >c d,ef [28] 

1 3 

W2 — — L a b, a b,ab — 3L a b !a b,ac + —L a b,ab,cd J T 
-\SL a b,ac,bd ~\~ 2L a b : ac,ad 6-L a c, bc.de 4" 2£/ a b )C d,e/ • [29] 

From the expressions for w\ and W2 we can extract the gen- 
eral expression for the exponent parameter A. However all 
the calculation above rely on the assumption that the repli- 
cated Gibbs free energy can be computed exactly. This is not 
possible in the general case and, as we have said in the main 
text, we have to recast in some given mean field-like approx- 
imation which has the correct glassy phenomenology. Here 
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we will give all the expressions above in the HNC approxi- 
mation where the derivatives of the Gibbs free energy can be 
computed exactly. The expressions for a and p axe 



2k 

i 

(7 = 

Sp 



(2tt) d 
{2-k) d 



kl(q) [l-pAc(q)]-K 



k 2 (q) [pAc(q) - 1] x 



' p 2 g 2 (x) 



Ac"(q) 



Ad{q) 



2 Q , Ac'(g) 
cos H — 



s/^*8W(M«)) 



2 2 
COS 6 1 



[301 



where Ac(q) = c(g) — c(g) is the difference between the diago- 
nal and off-diagonal part of the matrix of the direct correlation 
functions defined through the Ornstein Zernike equation and 
6 is the polar angle in D-dimensional polar coordinates. The 
expressions for the other two mass terms is given by 



7712 = 



7713 = — 
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d D q 



[31] 



(27T) 



M<?)s (?) 



By computing the third derivative of the replicated Gibbs free 
energy in the HNC approximation we get the expression for 
«ii and W2- 
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Ginzburg Criterion. In this section we give a guideline for the 
computation of the Ginzburg Criterion. In the main text we 
have said that at the dynamical point where the number of 
replicas goes to one, the leading behavior of the correlation 
functions of the two points function p^ can be computed 
using a field theory for a scalar quantity described by a cu- 
bic potential in a random field. This observation simplify 
a lot the loop expansion because it does not involve replica 
indices that complicate the perturbative analysis. With ref- 
erence to the action defined in Eq. [18] of the main text, we 
can give a perturbative expression for the two point function 
of the field (j>(x). The bare propagator is given as usual by 
{p) = up + P\ft + 8m. To obtain the two point function 
it is quite useful to write down the generating functional of 
the correlation functions W[J] = lnZ[J] where we can put 
J(x) — ho(x) + 8h and ho is an external field that can be used 
to extract the correlation function by taking the derivative 
with respect to it. Introducing the following diagrammatic 
notation 



J{x) = o 
we have that 



ho(x) 



Sh(g,A) 



[331 
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We impose that the critical point is not shifted by the per- 
turbative terms so we want also that {4>{x)) = from which 
we see that the counterterm 8h is of order g. Now let us look 
at the one loop correction to the propagator. Using the fact 
that the expectation value of <j> is zero we obtain 



(4>(x)ct>(y)} = G (x-y)+. 



■ + . 



We are interested in the most infrared divergent diagrams (in 
the limit where T — > Td). This means that we can neglect the 
second diagram, and we can consider only the first one (this 
is exactly what happens in the perturbative expansion of the 
Random Field Ising Model). The inverse of the renormalized 
susceptibility reads 



= G~ 1 (p = 0) =ml + 8m- 
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where m = pyft. By taking the derivative with respect to 
mo we obtain 



dm R i i q A 9 2 
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By imposing that the second term on the right hand side is 
smaller than 1 and by computing the loop integral we get the 
expressions [19] and [20] of the main text. 



Details on the numerics. To produce the numerical values col- 
lected in the tables, we have solved numerically the HNC equa- 
tions in three dimensions. This is a quite easy task because 
such equations can be solved by an iterative Picard scheme. 
However the solution requires the use of Fourier transforms. 
Working in spherical coordinates thanks to the rotational in- 
variance of the system, we have two natural cutoffs. The first 
one fixes the maximal distance L (infrared cutoff), hence we 
only keep g(r) for < r < L. The other one is related to the 
precision with which we measure the position of the particles 
(utraviolet cutoff): the possible values of r are discretized in 
such a way that in the unit interval there are N equi-spaced 
possible positions so that the precision is 1/N. The data pre- 
sented in the tables is relative to the larger cutoffs that we 
have. In particular, the infrared cutoff is fixed to L = 16 
where the unit distance is the diameter of the particles or the 
interaction range of the potential. The ultraviolet cutoff is 
fixed at N = 256. A remark has to be done on the way we 
computed the critical point and the zero mode. In fact to 
observe the correct y/e behavior of the off-diagonal solution, 
we need to be quite close to the critical point because other- 
wise this behavior is hidden by the subleading e behavior. To 
give a precise estimate of the critical point we have collected 
a sequence of solutions of the HNC equation varying the tem- 
perature or the density, depending on the case under study, 
and we have fitted these data with a yfe behavior. Once we 
have identified the critical point we have computed the the 
zero mode using directly the definition given by Eq. [8] of the 
main text. 
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Table 1. Numerical values of the coefficients of the effective action and the physical quantities from the HNC approximation. 
For each potential, lengths are given in units of ro and energies in units of e, with fcs = 1. Data at fixed temperature, using 
density as a control parameter with e = pd — p- 
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0.0125 


SS-12 


1 


2.057 
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0.0000973 


-0.0116 


0.00132 


0.00378 
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14.2 
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-0.0418 


0.00397 
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Table 2. Same as table 1, but here the data are at fixed density, using temperature as a control parameter with 

e = T d - T. 
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1.27 
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0.00153 
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0.325 


0.00195 
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